Statistical mechanics for metabolic networks during steady-state growth

Which properties of metabolic networks can be derived solely from stoichiometric information about the network's constituent reactions? Predictive results have been obtained by Flux Balance Analysis (FBA), by postulating that cells set metabolic fluxes within the allowed stoichiometry so as to maximize their growth. Here, we generalize this framework to single cell level using maximum entropy models from statistical physics. We define and compute, for the core metabolism of Escherichia coli, a joint distribution over all fluxes that yields the experimentally observed growth rate. This solution, containing FBA as a limiting case, provides a better match to the measured fluxes in the wild type and several mutants. We find that E. coli metabolism is close to, but not at, the optimality assumed by FBA. Moreover, our model makes a wide range of predictions: (i) on flux variability, its regulation, and flux correlations across individual cells; (ii) on the relative importance of stoichiometric constraints vs. growth rate optimization; (iii) on quantitative scaling relations for singe-cell growth rate distributions. We validate these scaling predictions using data from individual bacterial cells grown in a microfluidic device at different sub-inhibitory antibiotic concentrations. Under mild dynamical assumptions, fluctuation-response relations further predict the autocorrelation timescale in growth data and growth rate adaptation times following an environmental perturbation.Comment: 12 pages, 4 figure

Adaptation to parasites and costs of parasite resistance in mutator and nonmutator bacteria

Parasitism creates selection for resistance mechanisms in host populations and is hypothesized to promote increased host evolvability. However, the influence of these traits on host evolution when parasites are no longer present is unclear. We used experimental evolution and whole-genome sequencing of Escherichia coli to determine the effects of past and present exposure to parasitic viruses (phages) on the spread of mutator alleles, resistance, and bacterial competitive fitness. We found that mutator alleles spread rapidly during adaptation to any of four different phage species, and this pattern was even more pronounced with multiple phages present simultaneously. However, hypermutability did not detectably accelerate adaptation in the absence of phages and recovery of fitness costs associated with resistance. Several lineages evolved phage resistance through elevated mucoidy, and during subsequent evolution in phage-free conditions they rapidly reverted to nonmucoid, phage-susceptible phenotypes. Genome sequencing revealed that this phenotypic reversion was achieved by additional genetic changes rather than by genotypic reversion of the initial resistance mutations. Insertion sequence (IS) elements played a key role in both the acquisition of resistance and adaptation in the absence of parasites; unlike single nucleotide polymorphisms, IS insertions were not more frequent in mutator lineages. Our results provide a genetic explanation for rapid reversion of mucoidy, a phenotype observed in other bacterial species including human pathogens. Moreover, this demonstrates that the types of genetic change underlying adaptation to fitness costs, and consequently the impact of evolvability mechanisms such as increased point-mutation rates, depend critically on the mechanism of resistance

Autoregulation of mazEF expression underlies growth heterogeneity in bacterial populations

The MazF toxin sequence-specifically cleaves single-stranded RNA upon various stressful conditions, and it is activated as a part of the mazEF toxin–antitoxin module in Escherichia coli. Although autoregulation of mazEF expression through the MazE antitoxin-dependent transcriptional repression has been biochemically characterized, less is known about post-transcriptional autoregulation, as well as how both of these autoregulatory features affect growth of single cells during conditions that promote MazF production. Here, we demonstrate post-transcriptional autoregulation of mazF expression dynamics by MazF cleaving its own transcript. Single-cell analyses of bacterial populations during ectopic MazF production indicated that two-level autoregulation of mazEF expression influences cell-to-cell growth rate heterogeneity. The increase in growth rate heterogeneity is governed by the MazE antitoxin, and tuned by the MazF-dependent mazF mRNA cleavage. Also, both autoregulatory features grant rapid exit from the stress caused by mazF overexpression. Time-lapse microscopy revealed that MazF-mediated cleavage of mazF mRNA leads to increased temporal variability in length of individual cells during ectopic mazF overexpression, as explained by a stochastic model indicating that mazEF mRNA cleavage underlies temporal fluctuations in MazF levels during stress

Time-lapse microscopy data

Nela Nikolic, Tobias Bergmiller, Alexandra Vandervelde, Tanino G. Albanese, Lendert Gelens, and Isabella Moll (2018) “Autoregulation of mazEF expression underlies growth heterogeneity in bacterial populations” Nucleic Acids Research, doi: 10.15479/AT:ISTA:74; microscopy experiments by Tobias Bergmiller; image and data analysis by Nela Nikolic

Data from: Increased gene dosage plays a predominant role in the initial stages of evolution of duplicate TEM-1 beta lactamase genes

Gene duplication is important in evolution, because it provides new raw material for evolutionary adaptations. Several existing hypotheses about the causes of duplicate retention and diversification differ in their emphasis on gene dosage, sub-functionalization, and neo-functionalization. Little experimental data exists on the relative importance of gene expression changes and changes in coding regions for the evolution of duplicate genes. Furthermore, we do not know how strongly the environment could affect this importance. To address these questions, we performed evolution experiments with the TEM-1 beta lactamase gene in E. coli to study the initial stages of duplicate gene evolution in the laboratory. We mimicked tandem duplication by inserting two copies of the TEM-1 gene on the same plasmid. We then subjected these copies to repeated cycles of mutagenesis and selection in various environments that contained antibiotics in different combinations and concentrations. Our experiments showed that gene dosage is the most important factor in the initial stages of duplicate gene evolution, and overshadows the importance of point mutations in the coding region

Data from: Increased gene dosage plays a predominant role in the initial stages of evolution of duplicate TEM-1 beta lactamase genes

Gene duplication is important in evolution, because it provides new raw material for evolutionary adaptations. Several existing hypotheses about the causes of duplicate retention and diversification differ in their emphasis on gene dosage, sub-functionalization, and neo-functionalization. Little experimental data exists on the relative importance of gene expression changes and changes in coding regions for the evolution of duplicate genes. Furthermore, we do not know how strongly the environment could affect this importance. To address these questions, we performed evolution experiments with the TEM-1 beta lactamase gene in E. coli to study the initial stages of duplicate gene evolution in the laboratory. We mimicked tandem duplication by inserting two copies of the TEM-1 gene on the same plasmid. We then subjected these copies to repeated cycles of mutagenesis and selection in various environments that contained antibiotics in different combinations and concentrations. Our experiments showed that gene dosage is the most important factor in the initial stages of duplicate gene evolution, and overshadows the importance of point mutations in the coding region

Data from: Altruism can evolve when relatedness is low: evidence from bacteria committing suicide upon phage infection

High relatedness among interacting individuals has generally been considered a precondition for the evolution of altruism. However, kin-selection theory also predicts the evolution of altruism when relatedness is low, as long as the cost of the altruistic act is minor compared to its benefit. Here, we demonstrate evidence for a low-cost altruistic act in bacteria. We investigated Escherichia coli responding to the attack of an obligately lytic phage by committing suicide in order to prevent parasite transmission to nearby relatives. We found that bacterial suicide provides large benefits to survivors at marginal costs to committers. The cost of suicide was low because infected cells are moribund, rapidly dying upon phage infection, such that no more opportunity for reproduction remains. As a consequence of its marginal cost, host suicide was selectively favoured even when relatedness between committers and survivors approached zero. Altogether, our findings demonstrate that low-cost suicide can evolve with ease, represents an effective host-defence strategy, and seems to be widespread among microbes. Moreover, low-cost suicide might also occur in higher organisms as exemplified by infected social insect workers leaving the colony to die in isolation

Lack of cations in flow cytometry buffers affect fluorescence signals by reducing membrane stability and viability of Escherichia coli strains

Buffers are essential for diluting bacterial cultures for flow cytometry analysis in order to study bacterial physiology and gene expression parameters based on fluorescence signals. Using a variety of constitutively expressed fluorescent proteins in Escherichia coli K-12 strain MG1655, we found strong artifactual changes in fluorescence levels after dilution into the commonly used flow cytometry buffer phosphate-buffered saline (PBS) and two other buffer solutions, Tris-HCl and M9 salts. These changes appeared very rapidly after dilution, and were linked to increased membrane permeability and loss in cell viability. We observed buffer-related effects in several different E. coli strains, K-12, C and W, but not E. coli B, which can be partially explained by differences in lipopolysaccharide (LPS) and outer membrane composition. Supplementing the buffers with divalent cations responsible for outer membrane stability, Mg2+ and Ca2+, preserved fluorescence signals, membrane integrity and viability of E. coli. Thus, stabilizing the bacterial outer membrane is essential for precise and unbiased measurements of fluorescence parameters using flow cytometry